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Following a suggestion of H. E. Rowe, the first-order system of coupled 
differential equations urith the independent variables representing amplitudes 
of coupled modes (the so-called coupled mode equations) has been rigorously 
transformed into similar equations with the independent variables repre- 
senting the power in the coupled modes for the case where the amplitude 
coupling coefficient is a Gaussian random process with a white noise spec- 
trum and zero mean value. Coupled mode equations in power have proven 
useful, and this work provides formal justification and identification of the 
variables and parameters. 

The eqiiations may allow one to approach coupled-mode problems directly 
for statistical results. 

T. INTRODUCTION 

Consider the coupled line equations: 

h'(z) = -\\h(z) +jc(z)h{z) (1) 

/,'(«) =jc(z)h(z) - YMz). (2) 

These equations are useful in describing effects of coupling between 
a signal mode, represented by a complex wave amplitude Io(z), and a 
single spurious mode, represented by I\{z), caused by geometric imper- 
fections in a multimode transmission line. These equations may be de- 
rived in two ways from basic principles: ' ' ' direct conversion of Max- 
well's equations to generalized telegraphist's equations, or allowing 
discrete converters to become continuous. Exact solutions are known in 
only a few special cases, one of which is the case of constant c(z). Mak- 
ing use of this solution and assuming c{z) is a stationary Gaussian ran- 
dom function with a white noise spectrum S and zero mean value, it is 
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possible to derive coupled differential equations with the expected values 
of the power in I (z) and U(z) as independent parameters. That is 

P '(z) = (-2«d - S )P (z) + 80PM (3) 

P/(«) = SoPo(z) + (-2 ttl - S )Pi(2) (4) 

where 

PoGO = </.(«)/■*(!)>. 

1 1 1 

These equations are similar to equations used by S. E. Miller. 6 

II. DESCRIPTION OF MATHEMATICAL MODEL AND COMPUTATION 

Consider (1) and (2) and assume the h(z) are known and c(z) has 

i 

the constant value c in (z,z + Az) ; then I (z + As) are given by 4 

i 



7o(z + Az)' 



= (e _r ° A, )G 



_J,(2 + Az) 

where G is a two-by-two matrix with elements 
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+ (^ a )«p[f (i-i«)^ 



^ 2 = ( a -^)exp[f d+,a)A,] 






— (1 - ja)A« 



?12 = ff21 = (^AT ,lxp 



AT 

"9" 



(1 + ja)Azl 



-(^r) exp K (1_ia)A2 . 



where 



a = 



Ar = r - Ti = Aa + jA/3. 
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Now compute (I oh*) and (IJ\*) from 



/o/o 

US 

/i/o" 



- ( e - 2QoAj )(G xg*; 



_ J+iJ 



[/o/o*l 


W 


/i/o* 


_/,/,*_ 



where 



where X denotes the Kronecker product (direct product) of the matrices 
G and G*. 6 If we compute ensemble averages, assuming the random 
function c(z) is independent of previous values and (c(z)) = 0, we have 

P (z + ta) = e- 2a °* ! [AP (z) + BP 1 (z)] (5) 

P x {z + Az) = e- 2a °*'[BP (z) + DP x {z)\ (6) 

Po(z) = (h(z)h*(z)) 
i i i 

A = (01,011*) 

B = (012012*) 
D = (022022*) 

since (011012*) - (0i20n*) = (012022*) = (022012*) = 0, being odd functions 
of c. Expanding in powers (c ) and Az and requiring 

lim (c )Az = £0 

<c»>-.» 

Az-.0 

and after considerable algebraic manipulation, which we omit, we obtain 
(3) and (4) 

P (z) = (-2ao - S )P (z) + SoPxiz). 
1 110 

These equations may be related to Miller's equations by replacing our 
Po, Pi, 2a , 2ai , and S by his Pi + P„ , P x , a lh + a nh , a xh , and a x% . 
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